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We calculate the dispersion relation associated with a solitary wave in a quasi-one-dimensional 
Bose-Einstein condensate of atoms confined in a harmonic, cylindrical trap in the limit of weak and 
strong interactions. In both cases, the dispersion relation is linear for long wavelength excitations 
and terminates at the point where the group velocity vanishes. We also calculate the dispersion 
relation of sound waves in both limits of weak and strong coupling. 
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Bose-Einstein condensates of trapped alkali-metal 
atoms offer a rich source of interesting non-linear phe- 
nomena. At mean-field level, the effects of atom-atom 
interactions can be described as a one-body potential 
proportional to the local density of atoms. The order 
parameter, i.e., the condensate wavefunction, then satis- 
fies a non-linear Schrodinger equation which also includes 
the effect of the confining (harmonic) potential. 

Many authors have discussed the properties of solitary 
waves in a Bose-Einstein condensate of trapped alkali- 
metal atoms [^]-[)) . Solitary waves were created and stud- 
ied experimentally by Burger et al. JTo[ ] and by Denschlag 
et al. pd| . Solitary waves result as a balance between 
the energy cost associated with the Heisenberg princi- 
ple, fi, 2 /2m£ 2 , where m is the atom mass and £ is the 
characteristic length of the solitary wave, and the energy 
gain due to the local density variation of the system, 
which is of order nUo . Here n is the atomic density, and 
Uq = 4:Tth 2 a/m is the effective two-body interaction ma- 
trix element, with a being the scattering length for atom- 
atom collisions. In the present study, we consider repul- 
sive interactions with a > 0. Equating these two terms, 
one sees that the characteristic size of a solitary wave 
is set by the coherence length £ which satisfies the equa- 
tion h /2m£ 2 = nUo. In the experiments of Refs. jujfy], 
the effective interaction between the atoms was repulsive, 
and the solitary waves were thus density depressions. For 
an attractive effective interaction, the solitary waves are 
expected to be elevations in the density. 

Since the atoms are confined, momentum is not a good 
quantum number. However, it is possible to use cigar- 
shaped traps which are very long along the z axis. These 
systems are quasi one-dimensional p|,p!2|, and the mo- 
mentum along this axis is approximately a good quantum 
number. 

An interesting question arises in this context. Some 
40 years ago, Lieb considered a purely one-dimensional 
Bose gas of atoms interacting via a contact potential and 
predicted two distinct modes of excitation |L3| . One was 
identified as the usual Bogoliubov mode. The other class 
of excitations was later shown by Kulish et al. Ifi j J to be 
associated with solitary waves (see also Ref. |15||.) These 
authors demonstrated that the dispersion relation result- 



ing from solitary wave excitation is associated to that 
predicted by Lieb. It is thus reasonable to ask (at least 
in the case of quasi one-dimensional atomic condensates) 
whether this "Lieb mode" exists. Actually, in the recent 
study of Ref. [jl6[ the Lieb mode was examined in one 
dimension. Although the theoretical prediction for this 
mode seems firm, it has never been observed experimen- 
tally. In this regard, it is interesting that Stamper-Kurn 
et al. and Ozeri et al. |D| have recently managed to 
probe the long wavelength phonon spectrum associated 
with the Bogoliubov mode in a cigar-shaped condensate 
of atoms using Bragg spectroscopy. As we argue below, 
the Lieb mode should be present in such a system and 
may be observable. Since for long-wavelength excitations 
the Lieb mode coincides with the usual Bogoliubov mode 
of sound waves, it is crucial that the momentum imparted 
to the cloud be appropriately large for the two modes to 
have distinct energies. One, for example, could excite the 
cloud using the method of phase imprinting, in order to 
create a solitary wave, and then measure the excitation 
energy and the corresponding momentum. 

In the present study we derive the dispersion relation 
associated with the Lieb mode. Reference (and, 
recently, Ref. pfj[ ) has dealt with the same problem 
for a different range of parameters using a full three- 
dimensional numerical calculation based on the nonlin- 
ear Gross-Pitaevskii equation. In this calculation, the 
solitary wave is found to be a hybrid between a one- 
dimensional soliton and a three-dimensional vortex ring, 
but the method is applicable in our limit, as well. Since 
the present results deal with weaker interactions, we have 
chosen to adopt a description which neglects the con- 
tribution of vortex rings. We distinguish between two 
limits. In the limit of weak interactions, hqIIq <C htoj_, 
where uq is the maximum density of atoms far away from 
the wave, and co± is the frequency of the trapping poten- 
tial transversely to the long axis of the trap, the result- 
ing equation is the ordinary nonlinear Gross-Pitaevskii 
equation. In the opposite limit of strong interactions 
tiqUq huj±, the resulting equation is a modified Gross- 
Pitaevskii equation, in which the nonlinear term is pro- 
portional to the magnitude of the order parameter. Fi- 
nally, we calculate the usual Bogoliubov mode in both 



1 



regimes of weak and strong interactions and comment on 
the limits of validity of our study. 

Model. We start with a T — Bose-Einstein conden- 
sate of atoms confined in a cylindrical harmonic poten- 
tial, V = muj±(x 2 +y 2 )/2 and assume wave motion along 
the z axis. There is no confinement along the z axis, and 
away from the wave, there is a uniform density of atoms, 
n(x,y), which is independent of z. For a short-ranged 
atom-atom interaction, Vi n t(r — r') = UoS(r — r'), the 
Gross-Pitaevskii equation for the order parameter ^ has 
the form 



iUdtV = {-h 2 V 2 /2m+U \^\ 2 + V)V. 



(1) 



Following Ref. [Q, we assume that the transverse di- 
mension of the cloud is sufficiently small and the corre- 
sponding time scale sufficiently rapid that the transverse 
profile of the particle density can adjust to the equilib- 
rium form appropriate for the instantaneous number of 
particles per unit length. The problem becomes one- 
dimensional, and the solitary pulse can be described by 
a local velocity, v(z), and a local density of particles per 
unit length, a(z) Q, a(z) = J dxdy \^f(x, y, z)\ 2 . With 
this assumption, the wavefunction may be written in the 
form \£(r,i) = f(z,t)g(x,y,a) ||, where g is the equi- 
librium wavefunction for the transverse motion which de- 
pends on time implicitly through the time dependence of 
a. We choose g to be normalized so that J \g\ 2 dxdy = 1 
and thus, from the equations above, \f\ 2 = a. 

To proceed, we consider two opposite limits, namely 
the weak-coupling limit and the Thomas-Fermi regime. 
The transition between the two limits occurs for a^a ~ 
1/4 H , where uq is the value of a far away from the wave. 

Weak- coupling limit. We first consider the weak cou- 
pling regime. Although this has traditionally been an 
academic limit, it is now possible to create Bose-Einstein 
condensates in cigar-shaped traps pl| which realize the 
weak-interaction limit. In this case \g\ 2 has a Gaussian 



form, \g\ 2 = (7roi) _1 e _ ^/ 0j -) . As shown in Ref. 
satisfies the equation 



ihd t f = -(h 2 /2m)d 2 J + hu>x(l + Mf\ 2 )f- (2) 

We see from this equation that / cx e"™ 1 l 1+2 °' To ' f as 
\z\ — ► oo. Thus, we rewrite Eq. (^) using the variable 
w = f e ^±d+2<"?o)t to btain 

ihd t w = -(h 2 /2m)d 2 w + hujj_2a(\w\ 2 - a )w. (3) 

Equation (||) has the standard (quadratic) nonlinear term 
and leads to a speed of sound, c w , which satisfies the 
equation mc 2 w — 2huj±aoa H [see also Eq. (|l5|).] Since 
(To = nctTra 2 ^, we see that mc~ = no~Uo/2. 

Writing w — \foe 11 ^ and separating the real and imag- 
inary parts of Eq. (|^) , we obtain the two hydrodynamic 
equations 
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era - mu ) — > ( 4 ) 



u(l - o- /a). 



(5) 



Here, we have imposed the boundary condition v — > for 
a — > ctq . The solution of Eq. (0) is 



a(z)/a Q - 1 



cos' 



cosh (z cos9/Q ' 



(6) 



where 9 = arcsin(u/cu,) and £ = 2£(no) with £(no) equal 
to the coherence length for no = ao/^OjJ (i.e., £ = 
aj_/(2croa) 1 / 2 ). The wavefunction w can also be written 
as 



= y^[zsin0 + cos0tanh(;zcos0/£)]. 



(7) 



Energy and momentum of the solitary wave. In the 
limit of weak interactions, Eq. (0) implies that 
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- 2huj j_ao~oww* 
+h(jj xaa 2 ,) dz, (8) 



where the final term, which represents the energy of the 
background density of atoms, ensures convergence of the 
integral. Equation (@) can be written as 



£ = 



2m \ dz J 

+hujj_a(o~ — do) 2 ] dz. 



h 2 a (d$_ 
2m \ dz 



(9) 



Since v = %d z <j>/m, Eqs. (Q) and (^) allow us to reduce 
Eq. (||) to the simpler form 



£ = 2huj±a J (cr — Co) 2 dz. 



(10) 



Given the solitary wave profile of Eq. (§), Eq. (0) yields 
£ = (aV2/3)£ cos 3 0, (11) 



where £q = hcj±(aoa) 1 ^ 2 aoa±. 
Calculation of the momentum V , 

V = —ih [ w* — — dz — m [ a(z)v(z)dz, 
J dz J 



(12) 



requires a more careful description of the solitary-wave 
profile for large \z\. This is most easily accomplished by 
requiring that the solitonic wavefunction is periodic on a 
large interval [-L/2, +L/2]. The solution to this problem 
for all L can be expressed analytically in terms of Jacobi 
elliptic functions |22j. For our purposes, it is sufficient 
to note the behaviour of a and cj> in the limit of large 
L. For z - O(L ), a is still given by Eqs. (§) and ((f). 
For positive u, Eq. (||) indicates that d z cj> is negative and 
that a net phase will accumulate as we move from the 
center of the solitary wave at z = towards z = +L/2. 
For z ~ O(L), however, the periodic a approaches a 
constant value larger than ctq by an amount proportional 



2 



to X/L, In this region, d z <f> is positive, the total phase 
accumulated between z = and +L/2 is precisely 0, and 
the periodic boundary conditions are fulfilled. Evidently, 
this large- z behaviour makes a finite contribution to the 
momentum in the L — > oo limit which is readily evaluated 
from Eq. (|l^). To leading order, the momentum can be 
evaluated using Eqs. (g) and (Q) to obtain 

f 00 f°° (<7 — Cq) 2 

V = m\ {a — ctq)v(z) dz = mu / dz. (13) 



Using the order parameter of Eq. (Q), we find that 



V = V (ttu/\u\ -20 -sin 20), 



(14) 



where "Po = anfi. The momentum "P was also de- 
termined in Ref. p5| using the macroscopic relation 
V = fu~ 1 (d£/du)du. This yields a result identical to 
Eq. (|lj) . Equation (|l4|) implies that a maximum momen- 
tum of 'Pmax = 7rcro?i is obtained for u = 0. Note that, to 
leading order in 1/L, the energy of Eq. ( pi] ) is unaltered 
by the imposition of periodic boundary conditions. 
Combining Eqs. ( |ll|) and (14), we arrive at the disper- 



sion relation £ (V) for the Lieb mode in the limit of weak 
interactions. The solid line of Fig. 1(a) shows this result. 
For V — » 0, £ = c w V, in agreement with the usual Bo- 
goliubov dispersion relation discussed below. The Lieb 
mode terminates at V = T^max where it has an energy of 
£/£ = 4^2/3. 

Bogoliubov mode in the weak interaction limit. Equa- 
tion (||) implies that the Bogoliubov mode obeys the dis- 
persion relation 



£o Vo\ \VoJ 8 V a 



(15) 



Choosing ana = 0.1, aj_ = 1 /im, and a = 100 A, the 
coefficient inside the square root is 125. This number 
is relatively large because the characteristic wavevector 
corresponding to Vn, i-e., <7o, is much larger than the in- 
verse coherence length, ^ 1 {n ). Specifically, ^'(fij) ~ 
(8cr a) 1/2 /a±- Thus, a £{n ) ~ a±(a /a) 1 / 2 ~ 10. The 
dotted line in Fig. 1(a) shows the dispersion relation of 
Eq. (|l^) for the above numbers. 

The strong coupling limit. In the limit of strong inter- 
actions, uqUq ^> huj± 1 we use the Thomas- Fermi approx- 
imation for the transverse profile, |g| 2 = 2(ttR 2 l )~ 1 (1 — 
p 2 /i? 2 L ), with R±/a± = 2{aa) 1 / A @. From Ref. § we 
again see that w — f e ~ 2lljJ ±l< 7 oa)t sa tisfies the equation 

ihdtw = ~(h 2 /2m)d 2 w + 2hcu ± a 1/2 (\w\ - \w \)w. (16) 

The effective equation obeyed by w now involves a mod- 
ified non-linear term proportional to \w\. Equation Jl6| ) 
implies a sound speed of mc 2 = huj±(aoa) 1 ^ 2 (see also 
Eq. (|22j).1 Since uqIIo = 2hui±(aoa) 1 ^ 2 , this becomes 
mc 2 = n Uo/2 @E|. 



Using Eq. (|T^) , we again write to = \fae 1 ^ and separate 
real and imaginary parts. The velocity is still given by 
Eq. (||). In addition, Eq. ( |l6| ) implies that 

ft 2 / gyg \ 2 1/2 3/2 1/2 3/2^ 

2^(^J =3^/(3^/ -3<7 a + c ) 

_ TOu2 ^_^l!. (17) 

2(7 

This equation provides a relation between u and the min- 
imum value of (7, cr m in. For a given u, a min is given by 
the non-trivial root (a ^ an) of the right side of Eq. p7|) . 

Energy and momentum in the strong interaction limit. 
Using Eq. ( |l5| ) we find that the momentum is given by 



V/V = u/c s / (y - l) 2 /y dx, 



(18) 



where x = z(aga)^ 4 /a± and y = a /an. Here, a is the 
solution of Eq. (|l7|) . In addition, Eq. ( |l6|) gives an energy 
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1/2 



3/2 



dz. (19) 



The final term in the integral again guarantees its conver- 
gence and corresponds to the energy of the background 
density. 

Since v = ?id z (f)/m 1 Eq. (|l7) and the formula v = u(l — 
ao/a) allow us to write Eq. (19|) as 



£ = (4/3)?k^a 1/2 J (2(7 3 / 2 - 3(7(7 



1/2 



a^ 2 ) dz. 



Introducing the unit of energy £ — huj^(aoa) 1 ^ 4: ao a 



(20) 



£/£ = - I {2y*l 2 



3y + 1) dx. 



(21) 



We have solved Eq. jl7| ) numerically to obtain a(z) for 
various values of u. This numerical solution was used in 
Eqs. ( ]lgj ) and ( ^l| ) to determine V{u) and £(it), respec- 
tively. The solid line in Fig. 1(b) shows the resulting dis- 
persion relation for the Lieb mode in the limit of strong 
interactions. The x axis is measured in units of Vq, and 
the y axis is measured in units of £ . As in the case of 
weak interactions, the slope of the curve for small values 
of V is equal to c s , and also the curve terminates at ttVq 
with £ sa 1.5 £ . 

Bogoliubov mode in the limit of strong interactions. 
Equation (16) implies that the dispersion relation asso- 
ciated with the Bogoliubov mode is 



£ 



i 



(22) 



For ana = 1, a±_ = 1 fim, and a = 100 A, the coefficient 
inside the square root is sw 2500, which is 3> 1, since the 
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characteristic wavevector corresponding to Vq, i.e., ag, is 
3> than £ _1 (no): <7o£(no) ~ <7 a±/ '(aoa) 1 / 4 . The dotted 
line in Fig. 1(b) shows the dispersion relation of Eq. (p3) 
for this choice of parameters. Once again, the Lieb and 
the Bogoliubov modes have identical dispersion relations 
in the limit of long wavelength excitations, as expected. 

The methods adopted here are expected to be reli- 
able for weak interactions || for which the width of the 
solitary wave is much larger than the transverse width 
of the cloud. In this regime the transverse degrees of 
freedom are frozen out, and the behavior of the sys- 
tem is essentially one-dimensional. As the strength of 
the interaction is increased, it has been demonstrated 
in Ref. || that the dark solitary waves become unsta- 
ble. (See also Ref. §.) As shown in this reference, 
dark solitary waves in a cylindrical trap become unsta- 
ble for noUo/huj± > 2.4. Assuming that the system is 
in the Thomas-Fermi regime, ao/^irR^) = no/2, with 
R 2 ± = 4a^_(<7oa) 1 / 2 jl2|], and therefore the corresponding 
critical value of a^a is « 1.2 2 — 1.44. Thus, our approach 
should also provide a reasonable variational description 
of the dispersion relation in the strong interaction regime 
provided only that dark solitary waves are stable. In ad- 
dition, we note that the Bogoliubov dispersion relation 
obtained here is in analytic agreement with numerical so- 
lutions [ p3[ in the large and small momentum limits for 
both weak and strong interactions. Differences of only a 
few percent for intermediate momenta are found for aoa 
as large as 1. 

In summary, we have calculated the dispersion rela- 
tion E = £(V) of a sound wave and of a solitary wave 
in a quasi one-dimensional Bose-Einstein condensate of 
atoms confined in a harmonic, cylindrical trap in the lim- 
its of weak and strong interactions. For solitary waves, 
in both limits the spectrum has the same qualitative be- 
haviour: it is linear for long wavelength excitations and 
coincides with the Bogoliubov mode. For shorter wave- 
length excitations, it lies below the Bogoliubov mode, 
and it terminates at a maximum momentum for which 
the group velocity vanishes. 
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FIG. 1. (a) Dispersion relation £ — £{P) in the limit of 
weak interactions for a solitary wave (solid curve) and for the 
Bogoliubov mode (dotted line.) The energy is measured in 
units of £o and the momentum in units of Vo- (b) Same as (a) 
in the limit of strong interactions, with the energy measured 
in units of £ . The inset shows the same graph on a smaller 
scale, for long-wavelength excitations. 
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